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We outline a method for characterizing deviations from the properties of a Standard Model (SM) 
Higgs boson. We apply it to current data in order to characterize up to which degree the SM Higgs 
boson interpretation is consistent with experiment. We find that the SM Higgs boson is consistent 
with the current data set at the 82 % confidence level, based on data of excess events reported by CMS 
and ATLAS, which are interpreted to be related to the mass scale riih ^ 124 — 126 GeV, and on 
published CLg exclusion regions. We perform a global fit in terms of two parameters characterizing 
the deviation from the SM value in the gauge and fermion couplings of a Higgs boson. We find 
two minima in the global fit and identify observables that can remove this degeneracy. An update 
for Moriond 2012 data is included in the Appendix, which finds that the SM Higgs boson is now 
consistent with the current data set at only the 94 % confidence level (which corresponds to ~ 2 cr 
tension compared to the best fit point). 



I. INTRODUCTION 



The main goal of the LHC physics program is to unravel the origin of electroweak symmetry breaking, 
with the leading explanation being the Standard Model Higgs boson. In considering the SM Higgs hy- 
pothesis, or more generally the correct description of electroweak symmetry breaking (EWSB) as realized 
in nature, it is important to synthesize the experimental results which are determined to date, in terms of 
symmetries that must be (at least approximately) respected in a successful effective description of physics 
at the weak scale. 

Probes of the mechanism of electroweak symmetry breaking can be summarized in terms of tests of 
the properties of the W^, Z bosons through electroweak precision data observables (EWPD), probes of the 
low energy flavour changing signatures of weak scale physics, and direct searches. The results of the first 
two of these probes indicate that a viable theory of EWSB can reduce parameter tuning by adopting an 
approximate SU(2)c custodial symmetry to comply with EWPD |[T]^, and imposing approximate Minimal 
Flavour Violation (MFV) ||3-|8|, in order to be robustly consistent with the results of the flavour physics 
program. The SM Higgs sector with renormalizable couplings to the fermions is exactly consistent with 
MFV by definition, but it is only approximately SU(2)c symmetric. 
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Recently, ATLAS and CMS have reported significant exclusion limits for a Standard Model Higgs 
boson over a broad range of masses. The current exclusion regions reported by ATLAS ||9l are 112.7 < 
ruh < 115.5 GeV, 131 < ruh < 237 GeV and 251 < mt < 468 GeV at 95% confidence level (CL). While 
the Higgs mass region excluded by CMS lITOll at this CL is 127 < m/j < 600 GeV. There is also a suggestive 
combination of experimental hints in multiple final state signatures pointing in both experiments to a Higgs 
mass mil — 124 — 126 GeV. Although certainly not conclusive, the evidence for a light resonance with 
this mass scale and approximately SM Higgs-like properties has also been augmented by the observation 
of 7 events (compared to an expected background of 2 events) within one GeV of 124 GeV in the pp 
h-jj ~^ 1133 channel at CMS ifTTTl . which would correspond to a SM Higgs boson produced through vector 
boson fusion. We will use these experimental results in this paper to determine to what degree current data 
is selecting a SM Higgs boson or not. 

These experimental hints, along with the lack of any clear evidence of new states discovered to date 
at the LHC are suggestive that an effective theory of the EWSB sector, including a light scalar resonance 
and the approximate symmetries of SU(2)c as well as MFV, is currently an appropriate description of the 
data. We will examine recent results from the LHC in this framework and ascertain up to what degree 
experiment is selecting a SM Higgs doublet at present. We perform a broad analysis of this question in 
Sections II and III in terms of two parameters characterizing the deviation from the SM value in the gauge 
and fermion couplings of a Higgs boson. We find two nearly degenerate minima in the global fit. The 
effective Lagrangian we employ emerges naturally in composite Higgs scenarios. We comment on their 
interpretation in terms of current data. 

The LHC results should be interpreted in the context of the indirect evidence for the SM Higgs and its 
properties in EWPD. We examine the consistency of the results of our global fit to LHC Higgs-like data and 
EWPD in Section IV. 

Finally, in Section V, we discuss how future measurements can be presented in a manner that can more 
precisely and efficiently refine the understanding of such an effective theory with the aim of clarifying if the 
correct description of EWSB is the SM Higgs mechanism or not. In particular we point out the utility of 
ratios of best fit signal strengths which can be provided by the experimental collaborations and which have 
the ability to resolve the degeneracy in the two best fit regions we find when examining LHC Higgs-like 
data. This is due to the fact that, although an individual experimental signal can be faked by a Higgs- 
like degree of freedom, it is difficult for such a degree of freedom to simultaneously reproduce the Higgs 



' These values for exclusions and the data used in the body of the paper correspond to pre-Moriond 2012 values. The update of 
the fit for Moriond 2012 data is in Appendix B. 
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predictions in other channels with different dependencies on effective couplings. 

II. THE EFFECTIVE THEORY 

We consider an effective Lagrangian with a light scalar resonance, denoted as h in the following, that 
includes the Goldstone bosons associated with the breaking of SU(2) x U(1)y U(1)q and the SM field 
content. A minimal description of these degrees of freedom is given by an effective chiral Lagrangian with 
a nonlinear realization of the SU(2) x U(1)y symmetry. The Goldstone bosons eaten by the W^, Z bosons 
are denoted by vr" , where a = 1, 2, 3, and are grouped as 

= 6^'^'^^"/^ (1) 

with V = 246 GeV. The field transforms hnearly under SU(2)l x SU(2)r as LT.{x) R) 

where L, R indicate the transformation on the left and right under SU(2)l and SU(2)r, respectively, while 
SU(2)c is the diagonal subgroup of SU(2)l x SU(2)r, under which the scalar resonance h transforms as a 
singlet. A derivative expansion of such a theory is given by |[T2l - [T4l 



C = ^id,hf-Vih) + -TriD,^^D^i:) 



l + 2a- + 6^ + 63^ + 

V V V"^ 



~{uldl)T. l + c^ + C2^ + --- I ■''^ " I+/I.C., with 




V{h) = ^-i/.^ + |(^)/.3 + |(r^)/.^ + .... (2) 

Here we have adopted for later convenience the notation of composite models lfT2l - [l4ll . Minimal Flavour 
Violation dictates that the sole source of flavour violation are the Yukawa couplings and that /i is a singlet 
in flavour space and that its couplings to fermions (c, C2, ...) are flavour-universal. (Note that MFV is also 
compatible with a proportional to a combination of the SM Yukawa matricies.) We adopt this assumption, 
although we note that the LHC data is essentially only sensitive to flavour violation linked to the large top 
coupling (scaled by c from its SM value) at present. 

The SU(2) X U(1)y ^ U(1)q subgroup of SU(2)l x SU(2)r is weakly gauged in such theories. We 
will fit the current data including the leading hnear Higgs-like coupling effects in the Lagrangian of Eq. Q. 
We will neglect, however, dimension five operators h G^'^G^y, h W'^'^W^,^, h B'^'^B^ y and other higher 
dimension operators in the fit^, focusing on the effect of a and c. This can be justified by UV model 
building, for example in the composite Higgs case. In general, the coefficients a, 6, 63, c, C2, ... are arbitrary 



^ See Ref. 1151 for a discussion on the impact of tiiese operators on Higgs production. 
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numerical parameters subject to experimental constraints. Model building in the UV of this effective theory 
fixes relations between the parameters. 

The completion of the theory with a SM Higgs boson fixes a = 6 = c = (is = ^4 = 1 and 63 = C2 = 0, 
and higher order terms in the polynomial expansion of the h field vanish. In this case h becomes part of a 
linear multiplet 



reducing the theory to the SM Higgs boson Lagrangian where the growth in ^/s of the high energy scattering 
cross section Wl Wl — )• Wl Wl is exactly cancelled. For other choices of these parameters the growth 
in ^/s of high energy scattering processes is only moderated. First we will explore to what degree a light 
scalar resonance leading to the current experimental excesses reported by CMS and ATLAS is consistent 
with a SM Higgs boson interpretation or not.^ Note that the fermiophobic Higgs hypothesis is frequently 
interpreted to correspond to the singular point where c = 0, a = 1, while a more general interpretation of 
this hypothesis is c = with a unfixed, but constrained by current data. 



The experimental hints for the Higgs boson that have been reported to date in each channel individually 
provide only marginal evidence. The best fit values to the signal strengths ^ = {a x Br)/((T x Br)^^/ for 
particular Higgs masses are summarized in Table I. 

In Table I we also give in the second column the expected mass sensitivity of the various experimen- 
tal signatures as well as in the third column the best fit value for fi together with the la experimen- 
tal error and the 95% confidence level limit. If the mass resolution is not explicitly quoted, as in the 
pp — ^ ZZ* — > l^l^l^l^ and the pp — > WW* — ^ l^ul^v results given by ATLAS, we display as an 
estimate the mass sensitivity for these channels as quoted by CMS. Since the mass sensitivity is such that 
the best fit masses can reasonably be interpreted to overlap, we do not introduce a correction factor in the 
fit to shift to a common mass value. We also display, in column two, the local significance when given by 
the experiments. In column four we schematically list the leading sensitivity of the signal in terms of the 
cross section rescaling and the effective branching ratios to final states XX given by Brxx [a, c] . The actual 
production cross section rescalings are combinations including subdominant channels, that can introduce 
further dependence on (a,c), see Section III.A. We do not include a look elsewhere effect as we are per- 

^ Many proposals and discussions on establishing the correctness of the SM Higgs hypothesis exist in the literature, including 
Refs. fBitZSl . 




(3) 



III. FITTING ENHANCED CROSS SECTIONS 
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Channel [Exp] 


m/i [GeV] (Local Significance) 


M (Mi) 


Scaling to SM 


pj3^77 [ATLAS] 


126.5 ±0.7 (2.8cr) ES] 


2;° ? Ea (2.6) 


^ Br-y-y[a, c] 


pp^ ZZ* ^1+ i- 1+ r [ATLAS] 


126± - 2% (2.1 cr) lEg 


1.2ljj m (4.9) 


Br2z[a, c] 


pp^WW* i+ ui- V [ATLAS] 


126± - 20% (1.4 cr) ED 


1.2t°j ilTl (3.4) 


~ BrvKw[a, c] 


pp '-f'jjj [CMS] 


124 ±3% ifTolfTTl 




~ Br-y-y[a, c] 


pp -> 7 7[CMS, b, > 0.94] 


124 ±3% IfTolfTTl 




^ Br-y^[a, c] 


pp 77[CMS, b, i?;f'" < 0.94] 


124 ±3% IfTolfTTl 


2.iti:l[iij 


~' Br-y-y[a, c] 


pp 77[CMS, e, i?^"" > 0.94] 


124 ±3% IfTolfTTl 


0.0+2-9 111] 


'-^ Br..y^[a, c] 


pp-> 77[CMS,e,i?^"" < 0.94] 


124 ±3% HOKIT] 


4. lit? 1 111 


^ Br-y^[a, c] 


pp^ ZZ* £- £+ £- [CMS] 


126 ±2% (1.5 cr) ifTTIflSl 


o.5t;:° nni (2.7) 


~ Br22[a, c] 


pp^WW*^£+iye-i? [CMS] 


126±20%lfT0ll29ll 


o.7±°:^ nni (1.8) 


^ BrvyvK[a, c] 


pp^bb [CMS] 


124 ± 10% m 


i.2l}:^ Col (4.1) 


- Bi-^j,[a,c] 


pp-^Tf [CMS] 


124 ± 20%||10l 


0.8ii;7 nni (3.3) 


(? Br7.T=[a, c] 



TABLE I: Summary table of reported signatures with events related to the Higgs mass scale of interest (m/j ~ 124 
GeV) where excess events have been reported. 



forming a global fit for a particular mass value m/j f« 124 GeV. The CMS photon measurements are split 
up in accordance with Ref. 1 10]. The 77 events are classified by the conversion of the photon in the crystal 
- defining a parameter i^g - and their location in the detector, being endcap - e or barrel - b. The data we 
fit to also have associated exclusion curves. We take these exclusions into account by another procedure 
described in the text. Note that the r+r^ searches at ATLAS are included in the exclusion analysis but not 
fit to in the signal strength best fit as the corresponding experimental error is not available. 

The excesses of events with approximately the same mass scale in various channels are suggestive of a 
resonance, that could be interpreted as evidence of a light Higgs boson. We will assume that these excesses 
of events correspond to the same underlying physics and fit the data to discern the degree up to which the 
excesses are consistent with a SIVI Higgs boson interpretation.^ Our procedure to perform a global fit to the 
current data is as follows. 

First we fit to reported values of /u* including the deviations in the SM predictions by allowing the 
parameters a and c to deviate from their SM values of 1. We include in the fits the effects of modified 
production cross sections and branching ratios due to the rescaling of the SM couplings by the parameters 
a and c. In order to carry out these fits we are required to make a set of assumptions, which are summarized 
and discussed in the Appendix. We will illustrate the sensitivity of the fit to the various assumptions by 



The local significance quoted in Table I should be taken as a reminder to the reader of the very marginal situation of the statistical 
significance of these signals at present. 
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varying them in the results presented in Section III A 



For example, consider the case of the event yield used to construct each for pp ^ ^j, which in this 
discussion we will assume is only produced through gluon fusion. The SM prediction for the Higgs boson 
producing such events with an integrated Luminosity f Cdt can be schematically written as 

/ Cdt agg^hMh^ll) . (4) 

J cuts 

We will assume that the effects of the coefficients a,c are to simply rescale the number of events in the 
various signal cross sections. We neglect small shape differences in the differential distributions that could 
affect the event yield as (a, c) deviate from (1, 1). Then the integration over luminosity with phase space 
cuts will essentially cancel in the constructed theoretical prediction for the ratio /i*. The sole effect of the 
deviation from (1,1) for (a, c) will be to generate an excess/suppression of events compared to the SM 
prediction. This can be directly fit to the reported experimental best fit value for this ratio. This is the 
procedure we will adopt. 

We construct a measure for a two parameter fit in the following way. We define the matrix C as the 
covariance matrix of the observables, and A6i as a. vector of the difference in the observed and predicted 
value of the ratio, as a function of (a, c). The measure is then given by 

= (AOif {C-\, (AOj) . (5) 

The minimum Xmin determined, and the 65%, 90% and 99% best fit CL regions are given by Ax^ < 
2.1, 4.61, 9.21, respectively, for = Xmm + ^X^- The confidence level regions are defined by the 
cumulative distribution function for a two parameter fit. The matrix C is taken to be diagonal with the 
square of the 1 a theory and experimental errors added in quadrature for each observable in the diagonal 
element. As correlation coefficients are currently not supplied by the experimental collaborations, off- 
diagonal correlation coefficients are neglected. For the experimental errors we use the quoted 1 a errors 
on the reported signal strength. The errors 6^^ for the individual channels i are. made symmetric by taking 
6fi' = ^/{S^l'+)^ + (V-)2/^/2, where /X* = cr^/cTsjy.f- For theory predictions of the cross section values 
and related errors, we use the numbers given on the webpage of the LHC Higgs Cross Section Working 
Group [30] for ruh = 124 GeV and ^/s = 7 TeV.^ We symmetrize the total cross section error and propagate 
the error to get a theory error on /i. Taking 

J2i r'{a,c)a' x Bri(a,c) 
^ = — r^7"5 ' 



^ The theory errors are currently subdominant compared to the large experimental uncertainties. 
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with r*(a, c) the appropriate rescaling factor for each cross section and defining the error on each to be 
6a^, we determine the error 6fi by 



6n = fj. 



\ 



(7) 



For each search channel, ATLAS and CMS provide ifTOl [TTl l26l |29l an exclusion upper limit on 
each 'signal strength' so that ji^ > /x^ is excluded at 95% C.L. This value is reported in Table I. The final 
exclusion limit from combining all channels is also provided as a function of the Higgs mass, with a 
SM Higgs boson mass being excluded whenever < 1. We obtain such combined limits by a simple 
procedure, solving for /i i the equation 

? (S - A)^ " ? (A^l - A)' ^ ^ ' 

where /i is the average of the individual /i*, the measured signal strengths for each channel (which make 
larger when they are nonzero)^. When there are no excesses (or for the purpose of calculating expected 
limits) one simply sets /i* = 0, in which case our simple recipe combines in quadrature the limits from 
different channels. Again we neglect correlations in the measured limits on the signal strengths. In applying 
this procedure to our case, note that we have mapped the reported CLs exclusion curve as 

4m^/^'(«,c) (9) 

through rescaling the production cross sections and the decay widths. The exclusion regions are based on 
all released Higgs search analysis channels from the experiments given in December 2011 and updated in 
February 2012, and are comprised of the signals listed in Table I (plus the ATLAS h — r+r^ analysis). 
The results of the fit are dependent on the SM values used for the masses of the known particles, gauge 
coupling constants in the SM, etc. We summarize the SM inputs, which we have used, in Appendix A. 



A. Results 



Using the above described procedure we perform a global fit to the currently available data, with results 
shown in Fig.[T] The global fit results are combined with the exclusion contours in Fig. [2] The current data 
is consistent with a SM Higgs boson interpretation which lies on the 82% CL contour. The global quality 
of the fit is x^/d.o.f ~ 1. There are two best fit regions. For each signal strength to which we fit, there 
is a mapping from the SM point in (a, c) space, given by (1, 1), to a nearly degenerate point in terms of 



* We contrast this semi-empirical approximate formula with a more precise determination of the combined limit in Appendix B. 
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Fit to LHC Higgs like data, inclusive Fit to LHC Higgs like data, gg only 




FIG. 1: Global fit results in the (a, c) plane for all best fit a jusu values given by ATLAS and CMS. The SM Higgs 
boson with a linear realization of the EW chiral Lagrangian corresponds to the point (1,1) and is within the yellow 
90% CL region. The 65% and 99% CL regions correspond to the green and light gray shaded regions. The plot on 
the left includes all production channels, the plot on the right only includes the gg production channel in the signals 
with inclusive Higgs production, rescaled as described in the text. The jA.o.i ~ 0.99 for the left plot, while 
X^/d.o.f ~ 1.03 for the right plot. We also show two lines characterizing the relationship between the parameters a,c 
in the minimal composite Higgs scenarios of Refs. Il3ni32ll . red dashed (for MCHM4) and blue dashed (for MCHM5). 
In Eq. Q, rephasing of the Higgs field h ^ —h maps the best fit regions shown to other physically equivalent regions 
in the (a, c) space. 

Ax^ with differing values of (a, c). This mapping can be understood as following from the contours of 
constant production shown in Fig.[5]which show an approximate c -f-)- — c symmetry in the production of the 
Higgs and its subsequent branching ratios. This symmetry is broken by the interference of loop diagrams 
in various channels, with the strongest breaking of this approximate symmetry in the 7 7 final state. The 
best fit values for the two best fit regions in the inclusive plot are (0.99, —0.66, 1.99) and (1.2, 0.60, 2.6) in 
terms of (a, c, Xmin)- "^^^ point for comparison is (1, 1, 5.44). It is of interest to determine a means by 
which the best fit region degeneracy can be resolved. We discuss such an approach in Section V. 

Note that in Figs.[T|[2]we present results which we label as 'inclusive' or 'gg only'. These fits differ in 
how the production cross sections are treated. For 'gg only' in the signals with 77, r+ final states, we 
only rescale the dominant gluon fusion production channel, neglecting subdominant channels. While for 
results labeled as inclusive the following production channels are included for each signal: 

(i) W+ W~/ZZ production via gluon fusion and vector boson fusion, 
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Fit to LHC Higgs like data, inclusive 




a 



FIG. 2: Global fit results in the (a, c) plane for all best fit a/asM values given by ATLAS and CMS, taking into 
account all production channels. Also shown are the exclusion contours in the (a, c) plane determined by mapping 
the SM exclusion into an effective exclusion of [<j /asM)[o-^ c] < 1. Again the 65, 90, 99% CL regions coiTespond to 
the green, yellow and gray regions in the plots. The exclusion curve derived from ATLAS data is given by the red 
(dashed) line (with the region to the right of the line excluded), the CMS exclusion curve is the solid blue Une. 

(ii) 77 production via gluon fusion, vector boson fusion (VBF), tth production, and associated produc- 
tion with and Z. (Note that the 77 events that have associated jets, interpreted to come from 
VFB, are treated exclusively.) 

(iii) bb production is summed over associated h and h Z production, 

(iv) r+ r~ production is summed over gluon fusion, vector boson fusion, tth production, and associated 
production with and Z. Previous analyses at CMS are reported to only include the VBF initial 
state as the tagging jets are used to eliminate Drell-Yan Z — )■ r+ events. The updated results use 
a modelling of the Drell-Yan spectrum based on measurements of Drell-Yan produced Z — > /i+ ^~ 
events, so that the updated analysis is more inclusive. 

We do not find dramatic differences in the fits using the two different procedures, and thus are lead to 
consider our approximations used in performing the rescalings to be satisfactory. See Appendix A for a 
detailed discussion. 

As can be inferred from Fig. [T] and Fig.|2j the direct fit to the data and the exclusion curves are selecting 
the same region of parameter space. We also performed exclusive fits to the following subsets of data. We 
combined the 6 b and r+r^ data for a test of the fermion couplings and combined the W"*" W~, ZZ data 
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Fit to LHC Higgs like data, b.r Fit to LHC Higgs like data, W*W~, ZZ, VBF Fit to LHC Higgs like data, 77 




FIG. 3: Global fit results in the (a, c) plane for all reported best fit (j/(Jsm values given by ATLAS and CMS using 
subsets of the data. The conventions for the plots are the same as in the previous figures. Fitting the data in subsets 
degrades the fit. The left figure testing the bb and t+ data essentially is unconstrained with x^/d.o.f ~ 0, the 
middle plot with the combination of W+ W^, ZZ and vector boson fusion based 77 has ax^/do.f ~ 0.28 while the 
right plot based on all 77 excesses has a x^/d.o.f ~ 0.62. 

as well as the vector boson fusion based 77 data for a test of anomalous gauge couplings. Furthermore, 
we did a joint fit to the 77 best fit signal strengths. The results are shown in Fig. |3] As can be inferred 
from the figures, the data are insensitive to the fermion and gauge couplings in the rr and bb channels, 
the sensitivity in the combination of the W^, ZZ data as well as the vector boson fusion based 77 
data is more significant, with a -roughly- c symmetric region being selected and large enhancements of the 
gauge boson couplings being disfavoured. The results of the 77 fits are most interesting as the interplay of 
the top and gauge boson loops in the 77 decay amplitude select a skewed region in the (a, c) space. One 
can also see that the 77 data are most constraining at present. It cannot be overemphasized that the signal 
significances being fit to globally, and particularly when fitting to subsets of the data, are marginal. Strong 
conclusions are premature and the best fit values for x^/d.o.f in fitting to subsets of the data are poor. 

IV. INCORPORATING EWPD 

Deviations in the properties of a light scalar resonance from the properties of the SM Higgs boson also 
affect other low energy precision tests, in particular EWPD, which are sensitive to loop corrections to the 
gauge boson propagators due to Higgs boson loops as well as the SM gauge bosons loops. If the properties 
of a scalar deviate from the properties of the SM Higgs in the coupling to the SM gauge bosons, the high 
energy behavior of the theory is modified, and the decoupling of the longitudinal mode of the EW gauge 
bosons is no longer exact. 



11 



This effect manifests itself through the high energy behavior of the longitudinal degrees of freedom 
in high energy scattering still growing with energy, and in EWPD through loop diagrams involving the 
longitudinal degrees of freedom of the gauge bosons |33]. In EWPD, the corrections to the gauge boson 
propagators can be expressed in terms of shifts of the parameters STU Il34l436l given by 

Here we have introduced a cutoff scale A, which approximately represents the mass of the new states 
that are required in this framework to unitarize longitudinal gauge boson scattering at a scale given by 



A ~ 4 iTv/^/\l - a2|. For EWPD we use the results of the Gfitter collaboration |[37l 

5 = 0.02 ±0.11, r = 0.05 ±0.12, ?7 = 0.07 ± 0.12 . (11) 
And the correlation coefficient matrix is given by 

^ 1 0.879 -0.469^ 
C= 0.879 1 -0.716 . (12) 
y -0.469 -0.716 1 J 

We perform joint fits to LHC data and EWPD by adding the corresponding entries to an enlarged covariance 
matrix (including STU in terms of a) in the global fit. Correlations among EWPD observables are included. 
We then perform a new global minimization and joint fit. The results are given in Fig[4] We have included a 
correction to STU of the form [AS, AT, AU) = (0.003, -0.002, -0.0001) due to shifting the best fit value 
of the Higgs mass in these results from 120 GeV to 124 GeV using the exact expression for the one-loop 
Higgs boson contribution to STU. When a joint fit to current LHC data is combined with EWPD under 
the assumption that m/j ~ 124 GeV and with no other states in the EWPD fit, the tension in the fit slightly 
increases. The x^/d.o.f rises to 1.3, with EWPD favouring values of a > 1 as the STU parameters do not 
have vanishing central values for nih = 124 GeV. The allowed regions are consistent with the exclusion 
curves as shown in Fig [4] 



V. DISCUSSION 

In examining possible signal data for the Higgs boson, such as the reported best fit values of a/ asM 
from each experimental collaboration, the first task should be to clarify if a resonance is consistent with a 
SM Higgs interpretation or not. From an effective field theory perspective, there is no particular motivation 
for the couplings of the Higgs to be identical to one predicted by the SM. The hierarchy problem strongly 
suggests the existence of new states with a relatively low mass, coupled to the Higgs. Integrating out these 
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Fit to LHC Higgs like data + EWPD Data 



2- 
1 - 
c 0- 
-1 
-2- 

~(?.95 1.00 1.05 1.10 1.15 

a 

FIG. 4: Global fit results in the (a, c) plane for all best fit (j/(Jsm values given by ATLAS and CMS combined with 
EWPD. The plot follows the same convention as the previous plots. The x^/d.o.f = 1.3. As before the blue (red) 
curve is the exclusion region from CMS (ATLAS). Again the shown exclusion curves are determined from inclusive 
production. 

states, one expects the properties of the Higgs will change. Nevertheless, when approaching the data in 
effective field theory, the symmetries that are known to be at least approximately present at the weak scale 
are already highly constraining, as discussed in Section II. One is therefore lead to the effective description 
of the data which we have adopted^, and global fits to precision Higgs data will allow the SM hypothesis of 
EWSB to be experimentally tested in the near future. 

As the LHC will be running at 8 TeV, the contours of constant ^* = [cr x BR](/i — > XX)[a, c\/[a x 
BR]sAf(/i — ^ XX) in a final state X are constructed in the (a,c) plane for this centre of mass (cm.) 
energy in Fig. |5] We restrict ourselves to a subset of the channels used in this analysis. The prospects of 
reducing the degeneracy in the best fit regions relies mostly on further observations related to signal events 
in (J X BR(/i — 77). Here we have rescaled using the previously defined 'inclusive' rescaling and combined 
the production channels that contribute to the analyses of signal channels. Note that each individual channel 
has a mapping away from the SM point of (1,1) in the (a, c) plane to a family of degenerate points. This 
leads to two best fit areas, which have to be resolved. This can be done by realizing that the mapping is not 
the same when comparing among different channels. 

Experimental results can be presented in such a manner as to aid in this pursuit. When best fit signal 




^ The robustness of the requirement of exact MFV is less significant than SU(2)c in the EWSB sector. 
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Constant cr x 'Br{ZZ)/crsM ^ BrsM(ZZ) 



Constant cr x Br(yy)/o-sM x BrsM(rr) 
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FIG. 5: Contours of constant signal production in the (a, c) plane. We have also included the exclusion curves which 
we have derived from ATLAS and CMS data. 



strengths are presented for a common Higgs mass scale, correlation coefficients sliould be supplied as well, 
which will lead to more accurate fits. The effective theory approach makes clear that it is instructive to 
experimentally provide ratios of various best fit signal strengths, so that the parameters (a, c) can be more 
precisely determined. In particular, it would be helpful if the experimental collaborations provide results 
that allow the degeneracy in the two fit minima to be resolved. Ratios of effective signal strengths can help 
in resolving this degeneracy. In Fig.|6]we show the effective production contour ratios for two combinations 
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Fit to LHC Higgs like data, inclusive 
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Fit to LHC Higgs like data, inclusive 




^.0 




FIG. 6: Contours of constant ratios of signal production in the (a, c) plane overlapped on the best fit region. On 
the left, the ratio ji'^ / ^x^^ is superimposed for 8 TcV cm. energy on the best fit regions, on the right, the ratio 
I^vbfIP"''^ is superimposed. These observable ratios can resolve the degeneracy of the best fit regions. 



of extracted best fit signal strengths in the (a, c) plane. 



77 



(a X Sii)'^TT[a,c] (fTx5i?)|ff 



ISM 



I^Vef _ (o- X BR)^^p[a, c] (cr x BR) 



\h'yy 
SM 



,77 



{axBR)h^^[a,c] (a x BR)p^p 



(13) 



(14) 



SM 



Here we have included superscripts for the various production channels to make explicit the ratios to be con 
structed. For example a^^^ means the combination of the production channels discussed in Section 



III A 



that are included in ZZ signal events. These ratios are obviously 1 in the SM. It is important to note that 
comparing theoretical and experimental determinations of such ratios, which include sets of best fit sig- 
nal strengths simultaneously, will allow the degeneracy of the best fit regions to be significantly reduced. 
Such combinations can also be experimentally appealing when they allow systematic uncertainties to be 
cancelled, such as photon systematic uncertainties in p^y^p/ IjC"' ■ 



VI. CONCLUSIONS 



We have examined the current LHC data in an effective theory to determine to what degree the SM 
Higgs hypothesis is emerging from the data. To this end we have performed global fits of best fit signal 
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strengths and exclusion regions, taking into account current data. The SM Higgs hypothesis turns out to be 
consistent with the data at the 82 % CL. In our global fits we find that there are two best fit regions. We have 
determined experimentally accessible ratios of best fit signal strengths for a specific Higgs mass value that 
will allow the degeneracy in the best fit regions to be significantly reduced with sufficient data collected at 
8 TeV cm. energy. 

Appendix A: Fitting Assumptions 

In order to perform the fit we have made several assumptions, some of which can clearly be relaxed 
with more input from experimental collaborations. In this Appendix we discuss these assumptions in more 
detail. 

(i) We have assumed that the excess events reported in Table I, which have various best fit mass val- 
ues, correspond to the same underlying physics, which we are assuming can be characterized by the 
Lagrangian in Eq. Q. We have argued that symmetry considerations lead us to consider this La- 
grangian but the degree to which the reported best fit values of /i* can be associated with a particular 
mass scale is less clear. The mass resolution in all the various channels is larger than the spread 
of the Higgs masses listed in Table I, which differ at most by 2.5 GeV. However, reported best fit 
values of ^u* depend on the assumed rrih value. Considering the marginal statistical significance of 
the various channels, background fluctuations convoluted with signal events can be expected to shift 
the best fit value of the mass. Characterizing precisely the effect of such fluctuations is challenging. 
It is reasonable, however, to consider such effects to support our association of the excess with a 
common mass scale, which we choose to be 124 GeV. 

(ii) We have assumed that the given best fit values of ^u* in the various channels are uncorrelated to one 
another, neglecting both theory and experimental correlations. This is due to the lack of correlations 
reported by the experimental collaborations. This assumption is most easy to address with further 
experimental input. Our fitting procedure can be easily modified to include such correlation coef- 
ficients as off-diagonal elements in the covariance matrix. The signals naively expected to be most 
strongly correlated are the 77 events. We have studied the effects of correlations on the fit by intro- 
ducing pseudo-correlations through a correlation coefficient of size 0.5 between all 77 data. The fit 
results are robust against such pseudo-correlations and we still find two best fit regions with a similar 
parameter space as without correlations taken into account. We have also examined the robustness of 
the fit results against a set of other pseudo-correlations randomly chosen and assigned as off-diagonal 



16 



elements in the covariance matrix. Again we find robust fit results when globally fitting the data. 

(iii) We have assumed that the effect of rescaling the cross section and the branching ratios can be directly 
associated with a rescaling of neglecting the effect that rescaling the various channels modifies 
the differential distributions. This assumption can fail in the presence of higher dimensional oper- 
ators. The modification of the differential distributions will affect /x* due to experimental selection 
cuts modifying the shape of the differential distributions. However, we consider this effect to be 
subdominant to the effect we have incorporated by fitting to an unfixed (a, c) with production cross 
sections and branching ratios modified accordingly. Furthermore, we rescaled the various branching 
ratios according to the couplings included in the leading order formulae of the decay widths. This 
procedure is consistent when including only QCD corrections to the decay widths, however when in- 
cluding higher order EW corrections further effects due to a, c differing from 1 are neglected. Again, 
this effect is subdominant to the effects that we have retained. 

(iv) We have also neglected, in using /i*, that this result is determined assuming the SM in the combination 
of sub-channels, leading to the reported ^u* ratio. We also expect this effect to be subdominant to 
the effects that we have retained. This is in particular the case if the sub-channel combination is 
dominated by a particular final state. The analyses that use sub-channel combinations are the W W, 
ZZ, TT channels. We have tested the robustness of the fit by fitting with two procedures shown in 
Fig.[T] In the figure labeled 'gg only', we have only included the dominant production cross section 
leading to the various final states (usually (^gg^h production). We also have performed a fit where 
all production processed are included. The latter results are labeled as 'inclusive' in Fig.[T] We again 
find that the fit is robust against such differences. 

(v) We have neglected the effects of higher dimensional operators in the fit. This assumption can be 
justified in composite Higgs models, where the shift symmetry of the pseudo-Goldstone Higgs sup- 
presses such operators, and the coefficients we have retained are more sensitive to the underlying 
strong dynamics. 

(vi) The SM inputs used in the fit are nih = 124 GeV, mw = 80.398 GeV, mz = 91.1876 GeV, 
mt = 172.5 GeV and sin'^iOw) = 0.223. 

It is clearly of interest to further examine the effect of these assumptions, relaxing them whenever 
possible theoretically and to perform more refined analyses of this form when further experimental input is 
available. 
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Fit to LHC Higgs like data Fit to Global Higgs like data 




FIG. 7: Update to the global fit. To the left we have only updated the WW ATLAS signal in Table I to the value in 
Table II. On the right we have also added the ATLAS rr and Tevatron data on — > bb and pp ^ W^ W~ as shown 
in Table II. 



Channel [Exp] 


mh[GeV] 




VP^WW* ^ vl- V [ATLAS] 


126 


0.2t°:? (1.3) 


pp^bb [ATLAS] 


124 


-0.8^1;^ (3.5) 


pp^Tf [ATLAS] 


124 


(3.4) 


pp-^bb [CDF&D0] 


125 


2.0t!^;? (3-2) 


pp W+ W~ [CDF&D0] 


125 


0.03ti:22 (2.4) 



TABLE II: Summary table of reported signatures with events related to the Higgs mass scale of interest {mh — 
124 GeV) where excess events have been reported. Moriond 2012 update with new numbers to supplement (or 
replace) entiies in Table 1. 

Appendix B: Moriond 2012 Update 

In this section we present updated results including the data presented at Moriond 2012 1401 . The 
most significant change in the data that was used in Table I is an update to the ATLAS measurement of 
pp ^ W W* — )• £^ V £^ V. In addition, ATLAS reported best fit signal strengths in the pp ^ bb and 
pp ^ TT channels while CDF and D0 reported a broad excess in pp ^ bb events. Further, CMS has now 
also supplied best fit signal strengths as a function of rrih, allowing various Higgs mass hypotheses to be 
fit to. In this Appendix we include these experimental results in our fit and supply supplementary plots for 
various Higgs masses (refining also our determination of the 95% C.L. exclusion limits). 



18 



Channel [Exp] 


A119.5 (/^ 


119.5) 


A124 (/^ 


124) 


/il25 (M125) 


pp^ll [ATLAS] 


^•^-0.8 


(1.5) 


O.Strr 


(2.6) 


■■■•0-0.8 


(3.9) 


pp^ZZ*^e+ e- 1+ 1- [ATLAS] 


-0.5+1 


(5.1) 




(4.7) 


1 4+1-3 
■"■•^-0.8 


(4.1) 


VV^WW* ^ t^vl- V [ATLAS] 




(2.4) 




(1.6) 




(1.4) 


pp^^^ [CMS] 




(1.3) 




(3.5) 


1 6+°'^ 

■■^•0-0.6 


(3.0) 


pp^ZZ* ^(.+ 1- (.+ 1- [CMS] 




(5.2) 


0.5t;;^ 


(2.7) 


fi+°-^ 

U-D-0.6 


(2.5) 


W^WW* ^t^vt-v [CMS] 




(2.5) 


0.6t°;^ 


(1.8) 


0.41^6 


(1.5) 


pp^bb [CMS] 




(4.1) 




(5.0) 




(5.2) 


pp^Tf [CMS] 


o.2t?:? 


(3.6) 


0.4tl;0 


(3.9) 


o.6t\:l 


(4.1) 


pp^Tf [ATLAS] 


-0.9tl-l 


(2.9) 


-o.itl-J 


(3.4) 


n 1+1.7 


(3.5) 


pp^bb [CDF&D0] 


1 c+0.6 


(2.5) 




(3.1) 


^•"-0.7 


(3.2) 


pp W+ W- [CDF&D0] 




(4.5) 




(2.4) 




(2.4) 



TABLE III: Summary table of reported best fit signal strengths for various Higgs mass values. We note that the 
asymmetric nature of the error band for the best fit signal strength values reported for the pp ^ Z Z* —>■ £^ £~ £~ 
signal by ALTAS is curious. We use this data in the fit. For = 125 GcV we use the public results presented at 
Moriond 2012 that splits the 7 7 signal events into four classes (that are not identical to the classes used in the body 
of the paper) as well as the (dominantly) VBF induced photon events instead of the global photon value in the Table. 
The non VBF di -photon events we rescale inclusively. The data used for the category zero to four photons is given (in 
order) for ^125 by 2.1+1°^, 0.6^0.3, 2.21^;^ 0.5^1;^, while the VBF induced photon result is 3.6+1^. 

The updated data that we supplement Table I with is given in Table n. Due to an apparent inconsistency 
in the ATLAS best fit signal strength plot for pp ^ bb and the corresponding ATLAS CLg limit plot we do 
not use the b b best fit signal strength value in the combined fit. We show in Fig. 7 the effect of the Moriond 
2012 data on our previously reported fit results. 

We also show joint fits for the Higgs mass values m/j = 119.5, 124, 125 GcV where we have taken the 
experimentally reported fi and the corresponding theory predictions at a common m/^ due to the release 
of the required data by CMS, after version one of this paper. The data we use is given in Table III for 
each particular Higgs mass value chosen. These results are shown for various masses in Fig. 8. Examin- 
ing the results, one clearly sees that the excess at 119.5 GeV can be distinguished globally to be a likely 
statistical fluctuation compared to the global fit to the excess of events around ~ 124 — 126 GeV. For 
the three mass values, the best fit points {mh, a, c, x^) are (119.5, 0.28, 0.49, 9.8), (124, 0.87, -0.43, 4.1) 
and (125, 0.87, —0.42, 4.6), where for the latter two masses we use the CMS di-photon data resolved into 
subclasses. Fitting to rrih = 124 GeV with global di-photon best fit data only, compared to di-photon data 
split into subclasses, shows that the discrimination on the parameter space offered by separately reporting 
the VBF induced 77 jj signal is important. We also encourage the experimental collaborations to report the 
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mh =119.5 GeV mh =124 GeV, CMS yy split 




mh = 1 25 GeV, CMS yy split mh = 1 24 GeV, global CMS yy only 




l.O 0.5 1.0 1.5 2.0 l.O 0.5 1.0 1.5 2.0 

a a 

FIG. 8: Fitting various mass hypothesis using the data from Table III and its caption. The red dashed line is the ATLAS 
exclusion limit as described in the text, the blue solid line is the CMS limit and the combined CMS and ATLAS limit 
is included as a black solid line. 

degree of contamination of this signal with gg initial state Higgs events to enable a consistent treatment of 
the reported best fit signal stengths. One can demonstrate how the VBF signal interpolates between the two 
results shown in Fig. 8 by adding in contamination due to a{gg — )• h) events with our consistent rescaling 
procedure. One finds the series of plots shown in Fig. 9 for various degrees of contamination. Note that 
CMS also reports W~ jj events which offer a similar discrimination of the parameter space as the 
77 J j signal in principle. However, again contamination due to a{gg — )• h) events will exist and is not re- 
ported by the collaborations. We do not use this data as a separate channel at present. The update to the data 
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Fit to Higgs like data, VBF with gg 1% 
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FIG. 9: Various degrees of contamination of the VBF signal with gg ^ h events. 



has a small effect on the CL of the SM Higgs hypothesis compared to the best fit value of the current data. 
Assuming no contamination due to gg for VBF events one finds that our previously reported global fit with 
the Moriond 2012 data update (but without correcting to a single Higgs mass value in the experimental best 
fit signal strengths) has the SM hypothesis residing on a 94 % CL curve around the best fit value of (a, c). 
Assuming a 3% contamination of the VBF events due to gg, the SM hypothesis remains consistent with the 
data at 93 % CL compared to the best fit value. For a direct comparison, the Fermiophobic scenario with 
a = 1, c = is consistent with the data at 96 %CL for the same global fit. When a 3% contamination due 
to gg events for the VBF diphoton signal of CMS is assumed, the same Fermiophobic scenario is consistent 
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130 140 
mh (GeV) 




130 140 

(GeV) 



FIG. 10: 95% C.L. limits on ajasu for CMS (left) and ATLAS (right). The official curve is the black solid line. The 
red curve is obtained from the simple approximation of Eqn.([8]). The green curves are more precise determinations of 
the limit as explained in the text. 



with the data at the 88 %CL. We do not consider a Fermiophobic scenario to be favoured by the global data 
or the pattern of deviations from the SM in the current data set. With such marginal signal events, statistical 
fluctuations in the data are still present affecting the pattern of deviations. 

Finally, the 95% CL exclusion curves from ATLAS and CMS in the plots of this Appendix have been 
determined using a more precise method than Eq. ([8]l. In the same spirit of Ref. |[39]| . for each individual 
search channel i, we have first approximated the corresponding probability density function of the signal 
strength parameter /x by a Gaussian pi(/u) oc Exp[— (/i — /i*)^/ (2(7^)^^ J]. We have obtained the quantities 
/i' and o"obs,i trying to get the best approximation to the reported 95% CL in that channel, (obtained from 

..i 

the equation Jq ^ pi{fi)dn = 0.95). In the case of CMS channels, we use the approximation of Ref. ||391| 



(which uses (Jobs i ~ cr, 



cxp,'t 



A^L ea;p/l-96 and obtains /x by solving the equation that determines /u^. In 



the case of ATLAS data, we find better agreement with the reported limits by directly using /x* = /i* and 
Cobs,i as provided^. To illustrate the precision of our approximations to the exclusion limits we compare 



them with the official 95% CL limit on a/asM in Fig. 10 In the left (right) panel we show the CMS 



(ATLAS) limit, with the official curve in black. The red curve is the simple approximation of Eq. ^ and 
the green curves are more precise determinations of the limit as explained above. The dashed green line 
corresponds to an approximate determination of /x* and fiobs,* as in Ref. |[39l . The solid green line (only 



' This is not the case of the bb channel, where we find a discrepancy between the limit derived from /i' and (Tobs.i and the reported 
hmit. We therefore do not use this channel to extract the combined hmit in the (a, c) plane. Although this is a subdominant 
channel, the impact of such discrepancy on the SM exclusion around 115 — 120 GeV is not negligible. In the case of the Z ^ 41 
channel, for some ranges of Higgs masses the reported (negative) value of p, is hitting a boundary and we could not use it to 
reproduce well the experimental limit. In these cases we use the same approximation as for the CMS channels. 
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shown for ATLAS) uses = fi^ and o"obs,i as provided. 
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